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Abstract 

Future applications in environmental monitoring, delivery of services and transportation of goods moti- 
vate the study of deployment and partitioning tasks for groups of autonomous mobile agents. These tasks 
are achieved by recent coverage algorithms, based upon the classic methods by Lloyd. These algorithms how- 
ever rely upon critical requirements on the communication network: information is exchanged synchronously 
among all agents and long-range communication is sometimes required. This work proposes novel coverage 
algorithms that require only gossip communication, i.e., asynchronous, pairwise, and possibly unreliable 
communication. Which robot pair communicates at any given time may be selected deterministically or 
randomly. A key innovative idea is describing coverage algorithms for robot deployment and environment 
partitioning as dynamical systems on a space of partitions. In other words, we study the evolution of the 
regions assigned to each agent rather than the evolution of the agents' positions. The proposed gossip 
algorithms are shown to converge to centroidal Voronoi partitions under mild technical conditions. 

Our treatment features a broad variety of results in topology, analysis and geometry. First, we establish 
the compactness of a suitable space of partitions with respect to the symmetric difference metric. Second, 
with respect to this metric, we establish the continuity of various geometric maps, including the Voronoi 
diagram as a function of its generators, the location of a centroid as a function of a set, and the widely-known 
multicenter function studied in facility location problems. Third, we prove two convergence theorems for 
dynamical systems on metric spaces described by deterministic and stochastic switches. 

1 Introduction 

In the not too distant future, networks of coordinated autonomous robots will perform a broad range of en- 
vironmental monitoring and logistic tasks. Robotic camera networks will monitor airports and other public 
infrastructures. Teams of vehicles will perform surveillance, exploration and search and rescue operations. 
Groups of robots will enable novel logistic capacities in the transportation of goods and the delivery of services 
and resources to customers. New applications will be enabled by the ongoing decreases in size and cost and the 
increases in performance of sensors, actuators, communication devices and computing elements. 

In these future applications, load balancing algorithms will dictate how the workload is shared amongst and 
assigned to the individual robots. In other words, robotic resources will be assigned and deployed to competing 
requests in such a way as to optimize some performance metric. Remarkably, load balancing problems in robotic 
networks are often equivalent to robotic deployment and environment partitioning problems. For example, in 
surveillance applications, optimal sensor coverage is often achieved by partitioning the environment and assigning 
individual robotic sensors to individual regions of responsibility. Similarly, in the transportation of goods or 
delivery of services, minimizing the customer wait-time is equivalent to a multi- vehicle routing problem and, in 
turn, to computing optimal depot positions and regions of responsibility. 

•This work was supported in part NSF grant IIS-0904501, ARO MURI grant W911NF-05- 1-0219, and ONR grant N00014-07- 
1-0721 A preliminary and incomplete version of this work appeared in the Proceedings of the 2009 American Control Conference, 
pages 2228-2235, St. Louis, Missouri, USA. The authors gratefully acknowledge Prof. Andrew R. Teel for Remark 14.41 

t Francesco BuUo and Ruggero Carli are with the Center for Control, Dynamical Systems and Computation, University of 
California at Santa Barbara, Santa Barbara, CA 93106, USA, {bullo , carlirugjOengineering . ucsb . edu. 

■fPaolo Frasca was with the Dipartimcnto di Matcmatica, Politccnico di Torino, and is now with the D.I.I.M.A., Univcrsita di 
Salerno and with the Istituto per Ic Applicazioni del Calcolo-C.N.R., Roma, Italy, paolo.frasca@ginail.com. 



Motivated by these scenarios, this paper considers the two following interrelated problems. The deployment 
problem for a robotic network amounts to the design of coordination algorithms that lead the robots to be opti- 
mally placed in an environment of interest. Deployment performance is characterized by an appropriate network 
utility function that measures the deployment quality of a given configuration. The partitioning problem is the 
design of coordination algorithms that lead the robots to optimally partition the environment into subregions 
of interest; even here the objective is usually achieved through the design of appropriate utility functions. 

Literature review 

A wide range of literature is relevant to our investigation. In next four paragraphs we review previous work on 
centroidal Voronoi partitions, distributed robotic deployment, animal territory partitioning, and mathematical 
modeling tools. 

The "centering and partitioning" algorithm originally proposed by Lloyd [25j and elegantly reviewed in the 
survey [Tl] is a classic approach to facility location and environment partitioning problems. The Lloyd algorithm 
computes centroidal Voronoi partitions as optimal configurations of an important class of objective functions, 
called multicenter functions. Besides their intended application to quantization theory [18j . centroidal Voronoi 
partitions have widespread applications in numerous disciplines, including statistical pattern recognition |22| . 
mathematical imaging [15j , geometric optimization [3] and spatial resource allocation |12| . Recent mathematical 
interest has focused on convergence analysis |13| . bifurcation analysis of low dimensional problems [38], and 
anisotropic partitions [H], among other topics. 

Distributed and robotic versions of the Lloyd algorithm have been recently developed in the multiagent 
literature; see the survey [28] and the text [9l Chapter 5 and literature notes in Section 5.4]. We briefly 
review this growing literature in what follows. Generalized centroidal Voronoi partitions are shown in [T7] to 
be asymptotically optimal for estimation of stochastic spatial fields by sensor networks. In [lOj convergence 
results are obtained for simple basic interactions among mobile agents such as "move away from the closest 
other agent" or "move toward the furthest vertex of your region of responsibility." Convergence to centroidal 
Voronoi partitions is established in [4] for a class of communication-less sensor-based algorithms (related to the 
classic clustering work by MacQueen [27]). In [36j adaptive coverage controls are proposed for environments 
described by unknown density functions. Dynamic environments and corresponding dynamic coverage problems 
are treated in |21j . In [31j partitioning policies are shown to achieve optimal load balancing in vehicle routing 
problems, i.e., problems in which a robotic network provides service to customers that arrive in real time in the 
environment. 

Territory partitioning via competitive behaviors is a classic subject of study in behavioral ecology; see [T] 
for a comprehensive survey. For example, it is known [2] that the foraging behavior of conflicting colonies 
of red harvester ants (Pogonomyrmex barbatus) results in non-overlapping dominance regions that resemble 
Voronoi partitions. Non-overlapping dominance regions akin to centroidal Voronoi partitions are well docu- 
mented in [SI 1371 m] for the mouthbreeder fish (Tilapia mossambica). Territoriality behavior and competition 
among prides of African lions {Panthera lea) are discussed in [29]. Overall, numerous animal species achieve ter- 
ritory partitioning without a central coordinating entity and without synchronized communication, but rather 
relying upon asynchronous accidental interactions and stigmergy. To the best of our knowledge about biological 
and engineering multiagent systems, asynchronous territory partitioning has been barely studied, see [37] for 
introductory ideas about animal behavior, and mathematical models and analysis are lacking. 

To finalize the literature review, here is a synopsis of the mathematical ideas from diverse disciplines that 
we bring to bear on deployment and partitioning problems. First, we adopt the so-called gossip communication 
models in which only peer-to-peer asynchronous communication links are required. This communication model 
is widely studied in the wireless communication literature; example references include [231 [7]. Moreover, we 
consider control systems on a non-Euclidean state space; the interest for non-Euclidean spaces has a rich history 
in nonlinear control theory, dynamical systems and robotics, including the early work [8] and a recent application 
to multiagent systems [35] . Finally, we adopt various tools from topology and from the study of hyperspaces of 
sets [30]. 
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Statement of contributions 



This paper uncovers novel mathematical principles and tools of relevance in coordination problems and multi- 
agent systems. 

We tackle partitioning and coverage control algorithms in innovative ways. First, we design algorithms that 
require only gossip communication, i.e., asynchronous, pairwise, and possibly unreliable communication. Gossip 
communication is a simple, robust and effective protocol for noisy and uncertain wireless environments. Gossip 
communication may be implemented in wandering robots with short-range unreliable communication (see the 
illustrative motion coordination strategy in Appendix IX]) . Second, we propose a paradigm shift in coverage 
control and multicenter optimization. Classically [TOt \Tl\ [2T| [28l [3T| [36l [38] . the state space for the coverage 
algorithms are the agents' positions, i.e., as a function of the agents' positions the environment is divided into 
regions and regions are assigned to each agent. Instead, in our approach, the agents' positions are no longer a 
concern: the state space is the space of partitions of the environment and the algorithm dictates how to update 
the regions. 

Within the innovative context of gossip communication and partition-based mechanisms, we devise a novel 
algorithm for multicenter and coverage optimization. Our gossip coverage algorithm is a peer-to-peer version 
of Lloyd algorithm and aims to compute centroidal Voronoi partitions. Which robot pair communicates at 
any given time is the outcome of either a deterministic or a stochastic process. We also propose a modified 
version that restricts communication exchanges to adjacent regions and that has suitable continuity properties. 
Simulations illustrate that our algorithms successfully compute centroidal Voronoi partitions. 

To formally establish the convergence properties of our proposed gossip coverage algorithms, we perform a 
detailed mathematical analysis composed of three steps. First, we develop suitable versions of the Krasovskii- 
LaSalle invariance principle for dynamical systems on metric spaces described by deterministic and stochastic 
switches. Convergence to a set of fixed points is achieved under a uniform deterministic or stochastic persistency 
condition. Second, we establish the continuity of various geometric maps, including (1) the Voronoi diagram as 
a function of its generators, (2) the location of the generalized ccntroid as a function of a set, (3) the widely- 
known multicenter function studied in facility location, and (4) the gossip coverage algorithms. These continuity 
properties are established with respect to the symmetric distance metric in the space of partitions. Third and 
final, we study the topology of the space of partitions. With respect to the symmetric difference metric we 
prove the compactness of a relevant subset of partitions. Specifically, we focus on partitions whose component 
regions are the union of a bounded number of convex sets. 

In summary, relying upon our extensions of the invariance principle, the compactness of a subset of the set 
of partitions, and the continuity of the various relevant maps, we establish the convergence properties of the 
proposed gossip algorithms. In short, the algorithms converge to the set of centroidal Voronoi partitions under 
mild technical assumptions and under the assumption that the gossip communication exchanges satisfy either 
a deterministic or a stochastic persistency condition. 

Organization and notations 

The paper is structured as follows. In Section [2] we review multicenter optimization and coverage control ideas. 
In Section[3]we state our asynchronous territory partitioning problem, provide a solution via the gossip coverage 
algorithm, state the convergence properties of the algorithm and report some simulation results. The following 
Sections |4j [5] and [6] develop the mathematical machinery required to prove the convergence results. Section |4] 
contains the convergence theorems extending the Krasovskii-LaSalle invariance principle. Section [5] contains a 
discussion about the compactness properties of the space of partitions. Section [6] states the continuity properties 
of the relevant maps and functions and contains the proof of the main convergence results. Concluding remarks 
are given in Section [T] 

We let IR>o and R>o denote the set of positive and non-negative real numbers, respectively, and Z>o denote 
the set of non-negative integer numbers. Given A C M'', we let int(^). A, dA and diam(A) denote its interior, 
its closure, its boundary and its diameter, respectively. Given two sets X and Y , a set-valued map T : X ^ Y 
associates to an element of X a subset of Y. 
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2 A review of multicenter optimization and distributed coverage 
control 

In this section we review a variety of known results in geometric optimization and in robotic coordination. In 
Subsection 12.11 we review the notion of environment partitions and we introduce the multicenter function as 
a way to define optimal environment partitions and optimal robot or sensor positions in the environment. In 
Subsection 12.21 we review some distributed control algorithms for agent motion coordination and environment 
partitioning based on the classic Lloyd algorithm. 

2.1 Partitions, centroids and multicenter optimization 

We let Q denote an environment of interest to be apportioned. We assume Q is a compact convex subset of 
with non-empty interior. Partitions of Q are defined as follows. 

Definition 2.1 (Partition) An A^-partition of Q, denoted by v = is an ordered collection of N subsets 

of Q with the following properties: 

(i) U,e{i,...,iv} Vi = Q; 

(ii) int(ui) n int(vj) is empty for all i, j G {1, . . . , N} with i ^ j ; and 
(Hi) each set Vi, i G {1, . . . , iV}, is closed and has non-empty interior. 
The set of N -partitions of Q is denoted by Vn- 

Let p = {pi, . . . ,pn) G denote the position of N agents in the environment Q. Given a group of N agents 
and an A'^-partition, each agent is naturally in one-to-one correspondence with a component of the partition; 
specifically we refer to Vi as the dominance region of agent i e {1, . . . , A^}. 

On Q, we define a density function to be a bounded measurable positive function (j) : Q —> R>o and a 
performance function to be a locally Lipschitz, monotone increasing and convex function / : R>o —>■ R>o- With 
these notions, we define the multicenter function Hmuiticcntcr : Vn x — > R>o by 



This function is well-defined because closed sets are measurable. We aim to minimize TYmuiticenter with respect 
to both the partition v and the locations p. 



ticenter function Timuiticcntcr has numerous interpretations. Here we review two applications entailing robotic 
networks. First, in a vehicle routing and service delivery example given vehicles at locations pi, assume 
that f{\\Pi~q\\) is the cost incurred by agent i to travel to service an event taking place at point q. Events 
take place inside Q with likelihood (j). Accordingly, Timuiticcntcr quantifies the expected wait-time between event 
arrivals and agents servicing them. 

Second, in an environmental monitoring application \28^ . assume the robots aim to detect acoustic signals 
that originate and propagate isotropically in the environment. Because of noise and loss of resolution, the ability 
to detect a sound source originating at a point q from a sensor at position pi is proportional to the signal-to-noise 
ratio (which degrades with \\q — pi\\). If the performance function f equals minus the signal-to-noise ratio, then 
'^multicenter quantifies the expected signal-to-noise ratio and detection capacity for acoustic signals generated at 
random locations. □ 

Among all possible ways of partitioning a subset of W^, one is worth of special attention. Define the set of 
partly coincident locations Sn = {p € \ Pi = Pj for some i,j G {1, . . . ,N}, i ^ j}. Given p G \ Sn, 
the Voronoi partition of Q generated by p, denoted by V{p), is the ordered collection of the Voronoi regions 




(1) 



Remarks 2.2 (Locational optimization) 




{^^(P))z=v defined by 



V^{p)^{qeQ\\\q-p^\\ < \\q-Pj\\ for aU 



(2) 
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In other words, the Voronoi partition is a map V : {Q^ \ Sn) Vn- The regions Vi{p), i G {!,..., N}, are 
convex and, if Q is a polytope, they are polytopes. Now, given two distinct points qi and q2 in M'', define the 
[qi] q2) -bisector half-space by 

i?bi«cctor(gi;g2) = {<zeM''| Ik -gill < \\q-q2\\}. (3) 

In other words, the set iybiscctorCgi! Q2) is the closed half-space containing qi whose boundary is the hyperplane 
bisecting the segment from qi to q2- Note that bisector subspaces satisfy i?bisector('7i; '72) 7^ ^fbiscctor('3'2; 9i) and 
that Voronoi partition of Q satisfies Vi{pi, . . . ,pn) = Q H ( Cij^i i?bisector(j'i;Pj)) ■ 

Each region equipped with a density function possesses a point with a special relationship with the multi- 
center function. Define the scalar 1- center function Tii by 

Wife^)= / f{\\p-<l\\)mdq, (4) 

J A 

where p is any point in Q and A is a compact subset of Q. Under the stated assumptions on the performance 
function /, the function p Ti.i{p; A) is strictly convex in p, for any set A with positive measure (we postpone 
the proof to Lemma 16. ip . Therefore, the function p 7ii{p; A) has a unique minimum in the compact and 
convex set Q. We define the generalized centroid of a compact set A <Z Q with positive measure by 

Cd(A) ^ argmin{Hi(p; A)\p^ Q}. (5) 

In what follows, it is convenient to drop the word "generalized," and to denote by Cd(t;) = (Cd(wi), . . . , Cd(wAr)) G 
the vector of regions centroids corresponding to a partition v G Vjv- 

Remark 2.3 (Quadratic and linear performance functions) // the performance function is f{x) ~ , 
then the global minimum of Tii is the centroid (also called the center of mass) of A, defined by 

Cd{A) = (^jj{q)dqy' j^qmdq. 

If the performance function is f{x) = x, then the global minimum ofTii is the median (also called the Fermat- 
Weber center) of A. See |3 Chapter 2] for more details. □ 

Voronoi partitions and centroids have useful optimality properties stated in the following proposition and 
illustrated in Fig. [TJ 

Proposition 2.4 (Properties of TYmuiticenter) For any partition v G Vat and any point set p G \ Sn , 

Hmulticcntor(V^(p),p) < Wmulticontor(w, _p) , (6) 
Hmulticcntcr(w, Cd(w)) < Hmulticentcr(w, _p) . (7) 

Furthermore, inequality (|6]) is strict if any entry ofV{p) differs from the corresponding entry of v by a set with 
positive measure, and inequality ([7]) is strict if Cd{v) differs from p. 

The statements in Proposition 12.41 originate in the early work by S. P. Lloyd [25j; modern treatments are 
given in [T3] and [HI Propositions 2.14 and 2.15]. Proposition 12.41 implies the following necessary condition: if 
a pair {v,p) with p ^ 5jv minimizes TYmuiticontorj then p = Cd(w) and v = V(j>) up to a set of measure zero. 
Accordingly, the partitions that minimize Timuiticcntcr have the following property. 

Definition 2.5 The partition v G Vat is centroidal Voronoi if it has distinct centroids, that is, Cd{vi) ^ Cd{vj) 
for all j i, and v — V{Cd{v)). 
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Figure 1: Illustration of Proposition l2.4l The left figure shows a sample 2-partition v and point set p in a uniform 
square environment. The value of the cost function Timuiticcntor at {v,p) is diminished by either replacing v with 
the Voronoi partition generated by p (see central figure), or replacing p with the centroids of v (see right figure). 




Figure 2: The Voronoi partition and the corresponding Delaunay graph 



2.2 Prom Lloyd algorithm to distributed coverage control 

Here we consider a group of robotic agents with motion, communication and computation capacities and we 
review a coverage control algorithm that determines the motion of each robot in a group and the associated 
partition in such a way as to minimize Timuiticontcr- In what follows, we restrict our attention to d = 2, that is, 
we assume Q C K^. 

To explain in what sense our algorithms arc distributed, we introduce a useful graph. The Delaunay graph [111 
|9] associated to the distinct positions p G \ Sn is the undirected graph with node set {pi}^i and with the 
following edges: {pi,Pj) is an edge if and only if Vi{p) r\Vj{p) is non-empty. In other words, two agents are 
neighbors if and only if their Voronoi regions intersect, see Fig. [2l 

The coverage algorithm we consider is a distributed version of the classic Lloyd algorithm |14| based on 
"centering and partitioning" for the computation of centroidal Voronoi partitions. The algorithm is distributed 
in the sense that each robot determines its region of responsibility and its motion plan based upon communication 
with only some neighbors. Specifically, communications among the robots takes place along the edges of the 
Delaunay graph. The distributed coverage algorithm is described as follows. At each discrete time instant 
t G Z>o, each agent i performs the following tasks: (i) it transmits its position and receives the positions of 
its neighbors in the Delaunay graph; (ii) it computes its Voronoi region with the information received; (iii) it 
moves to the centroid of its Voronoi region. In mathematical terms, for t G Z>o, 

pit+i) = cd{vipm. (8) 

A variation of the function Hmuiticontor is useful to analyze this algorithm. We define the positions-based 
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multicenter function Tivoronoi : \ Sn R>o by 

N 

Hvoronoi(p) = Hmulticenter(l^(p),p) = / f (Wq - Pi\\)4>{q)dq- (9) 

Because of the compactness of Q, a continuity property, and the monotonicity properties in Proposition l2.41 one 
can show [21 Theorem 5.5] that Hvoronoi is monotonically non-increasing along the solutions of ([H]) and that all 
solutions of ^ converge asymptotically to the set of configurations that generate centroidal Voronoi partitions. 
Additional considerations about convergence are given in [13] . 

3 Gossip coverage control as a dynamical system on the space of 
partitions 

In this section we present the problem of interest, our novel gossip coverage algorithm and its convergence 
properties in Subsections 13.11 13.21 and 13.31 respectively. In order to reduce the communication requirements of 
our algorithm, we propose an adjacency-based and continuous algorithm in Subsection 13.41 Finally, we report 
some simulation results in Subsection 13.51 

3.1 Problem statement 

The distributed coverage law, based upon the Lloyd algorithm and described in the previous section, has some 
important limitations: it is applicable only to robotic networks with synchronized and reliable communication 
along all edges of the Delaunay graph (computed as a function of the robots' positions). In other words, the 
law ([HI) requires that there exists a predetermined common communication schedule for all robots and, at each 
communication round, each robot must simultaneously and reliably communicate its position. Note that the 
Delaunay graph, interpreted as a communication graph, has the following drawbacks: for worst-case robots' 
positions, a robot might have — 1 neighbors in the Delaunay graph and/or might have a neighbor that is 
arbitrarily far inside the environment. Therefore, each robot must be capable to communicate potentially to all 
other robots and/or to robots at large distances. 

Given this broad range of undesirable limitations, the aim of this paper is to reduce the communication 
requirements of distributed coverage algorithms, in terms of reliability, synchronization and topology. Here are 
some relevant questions that constitute our informal problem statement: 

Is it possible to optimize robots positions and environment partition with asynchronous, unreliable, 
and delayed communication? Specifically, what if the communication model is that of gossiping 
agents^ that is, a model in which only a pair of robots can communicate at any time? Since Voronoi 
partitions generated by gossiping and moving agents cannot be computed by gossiping agents, how 
do we update the environment partition? 

To answer these questions, the next subsections propose an innovative partition-based gossip approach, in which 
the robots' positions essentially play no role and where instead dominance regions are iterativcly updated. 
Designing coverage algorithms as dynamical systems on the space of partitions has the key advantage that one 
is not restricted to working only with Voronoi or anyway position-dependent partitions. 

Example 3.1 (The Lloyd algorithm in the partition-based approach) The distributed coverage algorithm ([5]) 
updates the robots' positions so as to incrementally minimize the function T^voronoi; while the environment par- 
tition is a function of the robots ' positions. In this paper we take a dual approach: we consider an algorithm 
that evolves partitions. From this partition-based viewpoint, the coverage algorithm is an iterated map on Vn 
and equation ^ is rewritten as v{t -I- 1) = V{CA{v{t))). □ 

3.2 The gossip coverage algorithm 

In this subsection we present a novel partition-based coverage algorithm in which, at each communication round, 
only two regions communicate. Recall the notion of bisector half-space from equation ([3|. 
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Gossip Coverage Algorithm 



For all t G Z>o, each agent i G {1,...,7V} maintains in memory a dominance region Vi{t). The collection 
{vi{0), . . . , vn{0)) is an arbitrary polygonal A^-partition of Q. At each t G Z>o a pair of communicating regions, 
say Vi{t) and Vj{t), is selected by a deterministic or stochastic process to be determined. Every agent k ^ 
sets Vk{t + 1) Vk{t). Agents i and j perform the following tasks: 

1: agent i transmits to agent j its dominance region Vi{t) and vice- versa 

2: both agents compute the centroids Cd{vi{t)) and Gd{vj{t)) 

3: if Cd{v^it)) = Cd{vj(t)) then 

4: Vi{t + 1) := v,{t) and + 1) :== Vj{t) 

5; else 

6: v,{t+l) := (i;,(t)Ui;,(t)) nffbisector(CdK(t));Cd(T;j(t))) 

iscctor 

In other words, when two agents with distinct centroids communicate, their dominance regions evolve as 
follows: the union of the two dominance regions is divided into two new dominance regions by the hyperplane 
bisecting the segment between the two centroids; see Fig. [31 As a consequence, if the centroids Cd(w,;(<)), 
Cd{v j{t)) are distinct, then {vi{t+ l),Vj{t + 1)} is the Voronoi partition of the set Vi{t) \Jvj{t) generated by the 
centroids Cd{vi{t)) and Cd{vj{t)). 




Figure 3: The gossip coverage algorithm. The left and right figure contain the initial partition and the partition 
after one application of the gossip coverage algorithm. In the middle figure wc show the two centroids and (with 
a dashed line) the segment determining the bisector half-space. 



We claim that the algorithm is well-posed in the sense that the sequence of collections {v{t)}t^z>o generated 
by the algorithm is an iV-partition at all times t. Indeed, it is immediate to see that the first two properties in 
Definition 12.11 are satisfied at all time if they are satisfied at initial time. Finally, at all times t, each component 
of v(t) is clearly closed and has non-empty interior because of the following geometric fact: there exists no 
half-plane containing the interior of a region and not containing the centroid of the same region. 

Now, for any i, j G {!,..., TV} with i ^ j, define the map Tij : Vn Vjv by 



where 



V, 



if Cd(?;,) = Cd{vj), 



I {vi, . . . ,Vi, . . . ,Vj, . . . , Vn), otherwise, 
Vi = {vtUvj) niJbisoctor(Cd(wi); Cd(uj)), ^^^^ 

Vj = {viUVj) nifbisoctor(Cd(Wj);Cd(u,;)). 

The dynamical system on the space of partitions is therefore described by, for t G Z>o, 

v{t + l)=T.Mt)), (11) 
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together with a rule describing what pair of regions is selected at each time. We also define the set-valued 
map T : Vn =^ Vjv by T{v) = {Tij{v) | i, j G {f , . . . , N}, i ^ j}. The map T describes one iteration of the gossip 
coverage algorithm; an evolution of the gossip coverage algorithm is one of the solutions to the non-deterministic 
set-valued dynamical system v{t -t- 1) £ T{v{t)). 

Remark 3.2 (Gossip disk-covering control) We believe that our gossip and partition-based algorithmic ap- 
proach is applicable to a broad range of coverage problems. For example, the worst-case multicenter function J 10^ 
is defined by Hworst(p, i') = niax,;g|i _jv} maxggt,. \\q —Pi\\- Maximizing Tiworst is equivalent to covering Q with 
N disks of smallest radius centered at p. As in Proposition \2.4\ for any v £ Vjv and p £ Q^\Sn, one can prove 
'Hworst {V{p),p) < 'Hv,orst{v,p) and Hworst{v , CC{v)) < Hwoj:st{v,p), whcrc CC{v) is the array of circumcenters 
of the components of v. Hence, a gossip coverage algorithm for Tiworst is designed by replacing centroid with 
circumcenter operations in (jlO[) . We leave this and further extensions to future works. □ 

3.3 Analysis results for the gossip coverage algorithm 

In this subsection we state the main analysis and convergence results for the gossip coverage algorithm. 

We begin by studying the fixed points of T and by introducing an appropriate cost function with monotonicity 
properties along T. Regarding the algorithm's fixed points, we extend Definition 12.51 as follows. A partition v 
is mixed centroidal Voronoi if, for all pairs {vi,Vj) with i ^ j, either CA{vi) = Cd(wj) or [vi,Vj) is a centroidal 
Voronoi partition of Vi Uwj, that is, Vi = [vi i^Vj) ni7biscctor(Cd(wi); Cd(wj)). 

Lemma 3.3 (Fixed points of T and centroidal Voronoi partitions) For i,j G {1, . . . , N}, j ^ i, denote 
the set of fixed points of Tij : Vn ^ Vn by Fij = {w G Vn \ Tij{v) ~ v}. The following statements hold: 

(i) rij^i Fij equals the set of mixed centroidal Voronoi partitions; and 

(ii) if V is a mixed centroidal Voronoi partition satisfying CA{vi) ^ GA[vj) for j ^ i, then v is centroidal 
Voronoi. 

Next, we define the partition-based multicenter function TYccntroid ■ ^ K>o by 

N 

'^centroid (^) — '^multicenter ('^; 

Cd{v)) = Y. f{\\q-Gd{v,)mq)dq. (12) 

Lemma 3.4 (Monotonicity of TYcentroid along T) For i, j G {1, . . . , N}, i ^ j, 

HccntroidCTy (w)) < 7icontroid(w), for all V G Vn, and 

HcentroidCTy (w)) < Hcontroid (w) , iff Tij [v) and V differ by a set of measure zero. 

The proofs of Lemmas l3.3l and l3. 41 consist of elementary manipulations and are omitted in the interest of brevity. 
In short, we have established that the function TYccntroid monotonically decreases along each Tij when away from 
fixed points, and that centroidal Voronoi partitions arc the fixed points of all Tij provided ccntroids are distinct. 

We now prepare to state the main convergence result for T . We need to introduce some useful properties 
for sequences of partitions and for switching signals. 

Definition 3.5 (Non-degeneracy) A sequence of N -partitions {w(i)}(gz>o is 

(i) (uniformly) distinct centroidal if there exists e > such that, for all t G Z>o and i,j G {1, . . . , N}, i ^ j, 

one has \\ Cd{v,{t)) - Cd(vj{t))\\ > e; 

(ii) (uniformly componentwise) non- vanishing if there exists e > such that, for allt G Z>o and i G {!,..., TV} , 
the Lebesgue measure of Vi(t) is greater than e; and 

(Hi) (uniformly componentwise) finitely convex if there exists G N such that, for all t G Z>o and i G 
{1, . . . , N], the set Vi(t) is the union of at most I convex sets. 



9 



Moreover, the sequence v is said to be non-degenerate if it is distinct centroidal, non-vanishing and finitely 
convex. 

For example, a sequence of partitions is finitely-convex if each component of each partition in the sequence 
is the union of a uniformly bounded number of polygons with a uniformly bounded number of vertices. 

Definition 3.6 (Uniform and random persistency) Let X he a finite set. 

(i) A map a : Z>o —>■ X is uniformly persistent if there exists a duration A G N such that, for each x £ X , 
there exists an increasing sequence of times {tA;}/cGZ>Q C Z>o satisfying tk+i — tfc < A and <7{tk) ~ x for 
all k e Z>o. 

(ii) A stochastic process a : Z>o X is randomly persistent if there exists a probability p S ]0,1[ and a 
duration A G N such that, for each x G X and for each t € Z>o, there exists k G {!,..., A} satisfying 

¥[a{t + k) ^ x\cT{t),...,a{l)] > p. 

We are now ready to state the main deterministic and stochastic convergence results for the gossip coverage 
algorithm. It is convenient to postpone to Section 16.21 the theorem proof and the definition of convergence in 
the space of partitions. 

Theorem 3.7 (Convergence under persistent gossip) Consider the gossip coverage algorithm T and the 
evolutions {w(t)}tGZ>o <^ defined by 

v{t + l)=T„(t~,{v{t)), forteZyo, 

where a : Z>o — > {(*, j) G {1, ■ • ■ i N}'^ \ i ^ j} is either a deterministic map or a stochastic process. Then the 
following statements hold: 

(i) if a is a uniformly persistent map, then each non- degenerate evolution v converges to the set of centroidal 
Voronoi partitions; and 

(ii) if a is a randomly persistent stochastic process, then each evolution v, conditioned upon being non- 
degenerate, converges almost surely to the set of centroidal Voronoi partitions. 

The statements of Theorem 13 . 71 relv upon the assumption of non-degenerate evolutions. It is our conjecture 
that, starting from generic polygonal partition, this assumption is typically satisfied. We will discuss some 
numerical evidence to this effect in the next section. 

Lemma 13.41 indicates how the function Ticcntioid plays the role of a Lyapunov function for the dynamical 
system defined by T. To provide a complete Lyapunov convergence proof of Theorem l3.7[ we set out to establish 
three sets of relevant results. First, we need to establish extensions of the Krasovskii-LaSalle invariance principle 
for set-valued dynamical systems over compact metric spaces. Second, we need to establish the compactness 
properties of the space of non-degenerate partitions. Third, we need to establish the continuity of the relevant 
geometric maps. These three topics are the subjects of Section |4l [5] and [6l respectively. 

3.4 Designing an adjacency-based and continuous algorithm 

The gossip coverage map T has the undesirable feature that it entails communication exchanges between any 
two regions. This communication requirement might be too onerous in some multiagent applications. Ideally we 
would like to require communications only between adjacent regions, that is, regions whose boundaries touch, 
or between "nearby" regions. We believe such communication requirements may be easily achieved in robotic 
networks and we detail a sample implementation for robots with range-dependent Poisson communication in 
Appendix [XI Additionally, we aim to design a continuous gossip coverage map. We require the modified map 
to be continuous for technical reasons: the invariance principles we adopt for the convergence analysis require 
continuity of the dynamical system. 
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Motivated by this discussion, we modify the map T to rely upon only adjacency-based communication and 
to be continuous. First, we introduce a pseudodistance notion between sets. Given closed A d Q and B C Q 
with non-empty interior, define 

pseudodist(A, B) = inf{||a - &|| | (a, &) £ int(A) x int(B)}. 

Second, we select a positive constant 6 <^ diam((5) and denote by : Vn Vn the modified gossip coverage 
map to be defined in what follows. For any i,j £ {1, . . . , N}, i ^ j, we give the following partial definition: 




Cd{vi) — Cd(t'j)l| = O) or ( pseudodist(vi, Uj) > 6), 
Cd{vi) — Cd{vj)\\ > S) and ( pseudodist(wi, Wj) = O). 



(13) 



Therefore, if either Cd{vi) and Cd{vj) coincide or the pseudodistance between Vi and Vj is larger than (5, then 
Tfj{v) = V, that is, the map Tfj leaves the partition unchanged. Additionally, if the pseudodistance between the 
regions Vi and vj is zero {vi and Vj are adjacent) and the distance between Cd(i',;) and Cd(vj) is larger than (5, 
thenT,'jiv) = T,,{v). 

Next, we consider partitions that do not satisfy cither of the two conditions in definition p3p . We define the 
unit saturation function sat ; M>o [0, 1] by sat(a;) = a; if .t G [0, 1], and sat(a;) = 1 if x > 1 and the scaling 
function /3y : Vat — > [0, 1] by 

Pij{v) = sat(|| Cd{vi) — Cd(uj)||/(5) ^1 — sat(pseudodist(ui,i;j)/(5)^. 

The first condition and the second condition in (fT3|l correspond precisely to Pij{v) = and Pij{v) = 1, respec- 
tively. For partitions v satisfying < (3ij{v) < 1, wc aim to define so as to continuously interpolate between 
the identity map and the map T; see Fig. U for an illustration. Let _R, = n iJbisector(Cd(wj), Cd(i',:)) and 
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Figure 4: Modified gossip between close but not adjacent regions (0 < f3i2{v) < 1). The bisecting line 7_l 
borders the set R2 = V2 n i?bisector(Cd(wi), Cd(w2)) that, in the map T, is assigned to vi (see Fig. [3]). According 
to instead, only the subset R2 C i?2 is assigned to vi. Loosely speaking, the "width" of R2 equals /3i2(w) 
times the "width" of Ra, where "width" of a set is the maximum distance from a point in the set to 7_l. 



Rj = Vj n ^^bisoctor(Cd(wi), Cd{vj)). Define the line 7_l = (?i/biscctor(Cd(w-,), Gd{v.i)) and 



Pi — a, point in int (Ri) that is maximally distant from 7^, 



Pj = a point in int (Rj) that is maximally distant from 



(14) 



Next, note that for each q £ RiU Rj there exists a unique line, say 7g, that is parallel to "fj_ and passes through 
q. Based on this notion, wc define 

Ri = {q e Ri \ dist{pi,jq) < (w) dist(pi, 7^) or dist(g,7_L) > dist(pj, 7^)}, 



R-j 



{q e Rj I dist(pj,7q) < /3y (w) dist(pj, 7^) or dist(g,7_L) > dist(pj, 7^)}. 



We can now complete the partial definition p3|l . For all v with < f3ij{v) < 1, that is, for all partitions not 
already dealt with in definition ()13|) . we define 



T^iv) = {vi, ...,{v,\ R,) U . . . , {vj \ Rj) U i?„ . . . , vn). 



it\i entry 



jth entry 
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As discussed for T, one can prove that the map : Vn =^ Vjv defined by T^{v) ~ {Tfj{v) \ i,j E {1, . . . , N}, i ^ 
j}, is well-posed and has the following properties. 

Theorem 3.8 (Convergence of modified gossip map) Consider the modified gossip coverage algorithm 
and the evolutions {v[t)}t£Z>a defined by 

v{t + l)=T'^^,^{v{t)), forteZ>o, 

where a : Z>o — > {(«, j) S {1, ■ • ■ , N}'^ I * 7^ j} ^•s either a deterministic map or a stochastic process. Then the 
following statements hold: 

(i) if (J is a uniformly persistent map, then each non-vanishing and finitely-convex evolution v converges to 
the set of mixed centroidal Voronoi partitions; and 

(a) if a is a randomly persistent stochastic process, then each evolution v, conditioned upon being non-vanishing 
and finitely convex, converges almost surely to the set of mixed centroidal Voronoi partitions. 

3.5 Simulation results and implementation remarks 

We have extensively simulated the partition-based gossip coverage algorithm T on a 2-dimensional polygonal 
environment with uniform density and performance function f{x) = x^. Simulations have been implemented 
as a Matlab program, using the General Polygon Clipper Library to perform operations on polygons. We 
adopted the following communication model: at each iteration, a region pair is chosen, uniformly at random, 
among all pairs of adjacent regions. Fig. [5] is an illustration of a typical evolution. 




Figure 5: An example simulation of the gossip coverage algorithm with uniform random edge selection. The 
environment Q is a rectangle with uniform density, centroids are computed with performance function f{x) = x^, 
and N = 6 regions compose the partition. Snapshots of an evolution are shown for t S {0, 20, 50, 100, 300}. One 
may verify numerically that the sequence converges asymptotically to a centroidal Voronoi partition. At t = 20 
one of the regions is disconnected. 

Our first numerical finding is that the gossip coverage algorithm appears to converge to centroidal Voronoi 
partitions from all initial conditions. This is the same property that the Lloyd synchronous coverage algo- 
rithm ([8]) is known to possess. In other words, our numerically-computed sequences of partitions always converge 
to centroidal Voronoi partitions ~ even though our theoretical analysis (1) requires a continuous interpolation 
from T to and (2) does not exclude convergence to degenerate partitions where some component regions 
might have coincident centroids, or might have empty interiors, or might be composed of "polygons with an 
infinite number of vertices," that is, arbitrary sets. 

A second numerical finding is that, throughout numerous sample executions, the resulting polygonal regions 
rarely have complicated shapes and large numbers of vertices. This is good news because of our assumption 
of finite convexity and because large numbers of vertices affect both the computation and the communication 
burden of the gossip coverage algorithm. 

Finally, it is possible, and we have observed it numerically, to have evolutions of the algorithm that, before 
converging to centroidal Voronoi partitions, have components with disconnected regions. We believe that there 
might be applications where it is desirable to maintain connectivity of the components of the partition and, 
therefore, we sketch here how to design a connectivity-preserving algorithm. Note that the update step of the 
partition-based coverage algorithm amounts to the exchange, among the agents, of a region, which consists 
in general of several connected components. In the connectivity-preserving algorithm, such components are 
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considered individually, and each of them is traded only if this can be done without loosing connectivity; if not, 
the component is kept by the previous owner. Numerical simulations indicate that such an algorithm leads to 
centroidal Voronoi partitions as well. 

4 On the Krasovskii— LaSalle invariance principle: set-valued maps 
on metric spaces 

In this section we consider discrete-time set-valued dynamical system defined on metric spaces. Our goal is to 
provide some extensions of the classical Krasovskii-LaSallc Invariance Principle; we refer the interested reader 
to [SlUnilMl for recent invariance principles for switched continuous-time and hybrid systems. 

We start by reviewing some preliminary notions including set-valued dynamical systems, continuity and 
invariance properties, and Lyapunov functions. On a metric space {X, d), where X is a set and d is a metric on 
X, a set-valued map T : X ^ X is non-empty if T{x) ^ for all x ^ X. An evolution of the dynamical system 
determined by a non-empty set- valued map T is a sequence {a;„}„gz-^(j C X with the property that 

Xn+\ e Tixn)-, n e Z>o. 

In other words, we regard a set-valued map as a nondeterministic discrete-time dynamical system. For set-valued 
maps we introduce notions of continuity and invariance as follows. A set-valued map T is closed at x X if, 
for all pairs of convergent sequences Xk ^ x and x'/^ —>■ x' such that x'^ G T{xk), one has that x' S T{x). 
Additionally, T is closed on C X if it is closed at all w G W . A set C X is weakly positively invariant for 
T if T{w) CiW IS non-empty for all w G W . A set W is strongly positively invariant for T if T{w) C W for all 
w (^W. 

We are ready now to state a Krasovskii-LaSalle invariance principle for set-valued maps defined on metric 
spaces. This result extends the Global Convergence Theorem in [26| to more general Lyapunov functions. Its 
proof follows the lines of the proof of Theorem 1.21 in [5], and is thus omitted. 

Lemma 4.1 (Krasovskii-LaSalle invariance principle for set-valued maps) Let{X,d) be a metric space 
and T : X ^ X be non-empty. Assume that: 

(i) there exists a compact set W Q X that is strongly positively invariant for T; 

(ii) there exists a function [/ : — > R such that U{w') < U{w), for all w and w' G T(w); 

(Hi) the function U is continuous on W and the map T is closed on W . 

Then there exists c G M such that each evolution of T with initial condition in W approaches a set of the form 

5 C\U^^{c), where S is the largest weakly positively invariant set contained in 

{w &W\3w' T{w) such that U{w') ^ U{w)]. 

In this paper, given the metric space (X, d), we deal with a set- valued map T : X ^ X defined by a collection 
of maps Ti, . . . , T„i : X ^ X via the equality T{x) = {Ti(x), . . . , Tm(x)} for x ^ X . For this kind of set-valued 
maps, closedness is related to continuity of ordinary maps. 

Lemma 4.2 Let Ti, . . . ,Tm : X ^ X be continuous on W C X. The set-valued map T : X ^ X defined by 
T{x) = {Ti{x), . . . , T„i{x)} is closed on W . 

Proof. Let Wn w and w'^^ u>' be a pair of convergent sequences in W, such that w'„ G T{wn). We claim 
that the continuity of Ti, . . . , T™ implies w' G T{w). 

Note that, by hypothesis, for all n G Z>o there exists i„ G {1, . . . , m} such that w'^ — Ti^ (wn). Because the 
set {1, . . . , m} is finite, there exists an index j G {1, . . . , m} that appears infinitely many times in the sequence 
{«Ti}nGN- Consider the subsequences {w„,} C {wn} and {w'„^} ^ {w'„}, such that w'„^ = Tj{w„i). Clearly, we 
have that w„^ — » vu and w'^^ w' , where from the continuity of Tj it follows that w' — Tj{w). Thus, w' G T{w) 
and the claim is proved. I 

The following result is a stronger version of Lcmma [4.H for a particular class of set- valued dynamical systems. 
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Theorem 4.3 (Uniformly persistent switches imply convergence) Let {X, d) be a metric space. Given 
a collection of maps Ti, . . . , Tm '■ X — > X , define the set-valued map T : X ^ X by T{x) = {Ti(a:), . . . , Tm{x)} 
and let {xn\n£Z>a '^'^ evolution of T . Assume that: 

(i) there exists a compact set W <^ X that is strongly positively invariant for T ; 

(a) there exists a function U : W —^M such that U{w') < U{w), for all w G W and w' € T{'w) \ {w}; 
(Hi) the maps Ti, for i £ {1, . . . ,m}, and U are continuous on W ; and 

(iv) for all i G {l,...,m}, there exists an increasing sequence of times {nk \ k G Z>o} such that Xn^+i = 
Ti{xrn.) and (rife+i — n^) is bounded. 

If Xq G W , there exists c G M such that the evolution {a;„}„gz^(, approaches the set 

{Fin---nF„-,)nu-\c), 

where Fi ~ {w G W \ Ti{w) = w} is the set of fixed points of Ti in W , i G {1, . . . , m}. 

Loosely speaking, (i) the compactness of a strongly positively invariant set, (ii) a monotonicity property for a 
Lyapunov function, (iii) continuity properties, and (iv) uniformly persistent switches among finitely many maps, 
together ensure convergence of each evolution to the intersection of the fixed points of the maps. 

Proof. [Proof of Theorem 14. 3| Let S be the largest weakly positively invariant set contained in 

{w &W\3w' G T{w) such that U{w') = U{w)} = FiU---UF^. 

Since T is closed by Lemma 14. 2[ the assumptions of Lemma 14.11 are met; hence there exists c G R such that 
U{xn) c and x„ ^ 5 n U^^{c). 

Let ^{xn) denote the w-hmit set of the sequence {xn \ n G Z>o}. If we show that ^{xn) C (Fi n • • • nF,„) n U^^{c), 
then the statement of the theorem is proved. We proceed by contradiction. To this aim, let x G S r\U^^{c) \ 
((Fi n • • • n F„i) n U^^{c)) and let {a;„,,}/igz>o be a subsequence such that Xn^ ^ ^■ 

Observe that for each x G 5* \ {Fi n • ■ • flFm), there exists a non-empty set C {1, . . . ,rn} such that, 
X = Ti (x) if i G Tg, and x ^ Ti (x) if i ^ Tg. By the continuity of the maps T^, there exists 6 G M>o such that, 
if i ^ Zg, then Ti{x) ^ x for all x G i?5(x) = {a; G | d(x, x) < (5}. Let now 

75 = min i min iu ix) - U (Tdx)) U > 0. 

iels |^2:GS4(£) V ' ' ' ' " / J 

By hypothesis, if i ^ Xj, then U{Ti{x)) < U{x) for all x G Bs{x). Hence, since Bs{x) is closed, and U and the 
maps Ti are continuous, we deduce that 7 > 0. 

Observe now that hypothesis (iv) implies the existence of a duration Z? G N such that every map T,;, 
i G {l,...,m}, is applied at least once within every interval [n,n + D[, for all n G Z>o. Consider the set 
{Tij-gj^; this is a collection of continuous maps having x as fixed point. Then, there exists a suitable e G K>o 
such that, given any r-uple (ji, . . . , jV) £ 21, r < D, we have that Tj^ o Tj^ o Tj^ o . . . o Tj^(x) G Bs{x) for all 
w G B^{x). 

Select now k such that the element .t„^ in the subsequence {xnh}hez>Q satisfies d {x„^ ,x)<e and U {x„^. ) — c < js- 

Let 

s = min{t G [1, D[ \3j ^ Is such that Xn^+t+i = Tj (a;„^+t)}. 

Observe that U{xn^+s) - c < and C/(a-„^+s) - U (Tj{xni^+s)) > is implying that U (T, {x„^+s)) < c. This is 
a contradiction. I 

Remark 4.4 An alternate proof of this theorem can be given by applying an invariance principle obtained in J3^ 
on an appropriately-designed dynamical extension of the discrete-time set-valued system. □ 

In Appendix|B]we show how to persistent switching assumption is necessary. Next, we provide a probabilistic 
version of the previous theorem. 
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Theorem 4.5 (Persistent random switches imply convergence) Let {X,d) be a metric space. Given a 
collection of maps Ti, . . . , T„i : X — > X , define the set-valued map T : X ^ X by T{x) ~ {Ti{x), . . . , Tm{x)}. 
Given a stochastic process a : Z>o — > {1, • . • , m}, consider an evolution {a;„}„gz>o ofT satisfying 

*^n+l — ^cj(n)('^n)- 

Assume that: 

(i) there exists a compact set W C X that is strongly positively invariant for T; 

(ii) there exists a function U : W —> M. such that U{w') < U{w), for all w W and w' G T(w) \ {w}; 
(Hi) the maps Ti, for i G {1, . . . ,m}, and U are continuous on W ; and 

(iv) there exists p G ]0, 1[ and fc G N such that, for all i G {1, . . . , m} and n G Z>0; there exists h E {1, . . . , k} 
such that 

¥[a{n + h) =i\ a{n), . . . , a{l)] > p. 
If xq G W, then there exists c G M such that almost surely the evolution {xn}nei,-^o cippi^oaches the set 

iFin---nFra)nu-\c), 

where Fi ~ {w G W \ Ti{w) = w} is the set of fixed points of Ti in W , i G {1, . . . , m}. 

Loosely speaking, (i) the compactness of a strongly positively invariant set, (ii) a monotonicity property for a 
Lyapunov function, (iii) continuity properties, and (iv) persistent random switches among finitely many maps, 
together ensure convergence of each evolution to the intersection of the fixed points of the maps. 

Proof. [Proof of Theorem 14. 5| If xq G W, then the stochastic process cr induces a stochastic process taking 
values in W. From now on, we restrict our attention to sequences {x,i},igz>o such that xo G W. In other words 
we assume that the sample space containing all the evolutions of our interest is given by 

A ^ {{x„}„6z>o \ Xn &W for aU n G Z>o}. 

Let S be the largest weakly positively invariant set contained in 

{w eW\3w' e T{w) such that U{w') = U{w)} = Fi\J-- • U-F,„. 

From Lemma 1331 we know that there exists c G K such that a-„ S'nJ7"^(c). This implies that the following 
event is certain: 

E = \ {xn}n£Zyn I 3c G M such that lim U{xn) ~ c|. 

Let fi(a;„) denote the w-limit set of the sequence {xn I n G Z>o}. If we show that ri(.T„) C ((Fi n • • • n i^„i) n [/^^(c)) 
almost surely, then the statement of the theorem is proved. Next, consider the event 

El = {{x„}„6z>o \3x e S\{Fin. . .DF^) such that x G n{x„) }. 

Assume by contradiction that P [£'i] > 0. Now we compute P [£'|i?i]. Note that, for each x G 5'\ (i^i fl ... n F™), 
there exists a non-empty set C {1, . . . , m} such that, x = Ti (x) if i G 2£, and x ^ Ti [x) if i ^ Xj. Similarly 
to the proof of Theorem 14. 3[ we can associate to each x a positive real number 5 such that the inequality 
X j^Ti {x) holds true for all x G Bs{x) = {x & W \ d{x, x) < 6} and for all i ^ Xj. Moreover, we can define 

7^ - min ( mill (u{x) - U{T,{x))) ] , 

where the continuity of the maps Tj, j G {!,... ,rn}, and [/, and the closedness of the set Bs{x) ensure that 
is > 0. 

Consider the set {TjI^^j^; this is a collection of continuous maps having x as fixed point. Therefore, there 
exists a suitable e G M>o such that, given any r-uple (ji, . . . ,jr) G T~, r < k, we have Tj-^ °Tj2 o. . .oTj^{x) G 
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Bs{x) for all x G B^{x). Given {x„}„gz>Q, if there exists x £ S* \ (Fi n . . . fl i^,„) such that x e (a^n), then 
there must exist {nh\ h G Z>o} such that G B^{x) for all /i G Z>o. Moreover, without loss of generality we 
can assume that rih+i — Uh > k for all h G Z>o. Consider now the event 

E3 = {{in e {!,..., to} I 71 G Z>o} I 

3 h such that in^+s G -^s for all s G {1, . . . , fc} and ft, > ft,}. 

To compute P[£'3], we define, for j G Z>o, 

-^3j = {{in e {1, . . . , to} I n G Z>o} I «ri,,+s G -^S for aU s G {1, . . . , fc} and h> j}. 

Observe that {E^j \ j G Z>o} is a countable collection of disjoint sets such that = Ujlo -^3.j- -l^y hypothesis 
we have that 

F[E3,,] < lim r[(l-p) = 0, 

and therefore P[-E'3] = 0. This implies that, almost surely, there exists a subsequence {nhs}sei.>o ^ {"/i}/igz>o 
with the property that, for all s G Z>o, Xn^^+i = Ti{xn^^) for some i ^ and, therefore, also with the 
property that U{xn^^) — C/(a;„^^+i) > 75. Consequently, almost surely, we have that lims_^oo U{xn^^) = —00 
thus violating the fact that is a certain event. This implies that P [Ei] = and that, almost surely, Xn 

(Fin...ni^™)n[/-i(c). ' ■ 



Remark 4.6 The assumption, in Lemma \4.1\ and Theorems \4-3\ and \4-.5[ that W is strongly positively invariant 
ensures that any evolution of T with initial condition in W remains in W . By relaxing this assumption, it 
is possible to provide weaker versions of these results. Specifically, requiring W to be only compact and not 
necessarily strongly positively invariant, the thesis of Lemma \4-l\ o,nd Theorems \4.3\ and |^.5| do not hold in 
general for any evolution of T with initial condition on W, but are still valid for those evolutions {x„},igz>o 
that take values in W for all n G Z>o. □ 



5 On the topology of the space of partitions: compactness properties 
in the symmetric difference metric 

Motivated by the invariance principles presented in Section|4l we study metric structures on the set of partitions, 
with a focus on compactness and continuity properties. Specifically, we show how a particular subset of the 
set of partitions can be regarded as a compact metric space and how certain relevant maps are continuous over 
that subspace. In this section, and only in this section, the assumptions on Q are relaxed to give more general 
results: we assume that Q C M'^ is compact and connected and has non-empty interior. 

Let C denote the set of closed subsets of Q. We would like to introduce a topology on C with two 
properties: C is compact and the Voronoi map, the centroid map, and the multicenter function, defined in 
equations ^ (O, and respectively, are continuous over C. A natural candidate is the topology in- 
duced by the well-known [35] Hausdorff metric on C: given two sets A,BgC, their Hausdorff distance is 
d}i[A,B) = max {maxog^ min^gs c?(a, 6), maxbgs min^g^ d(a, 6)}. This metric induces a topology on C which 
makes it a compact space, but is not suitable for our purpose because, with respect to this topology, the Voronoi 
map, the centroid map, and the multicenter hmction, are not continuous; see Appendix [Cl Additionally, note 
that, unlike the Hausdorff metric, the centroid map and the multicenter function arc insensitive to sets of 
measure zero. 

In what follows, we introduce the symmetric difference metric, as a metric insensitive to sets of measure zero. 
Given two subsets A,B €C, we define their symmetric difference by AqB = {AlJB)\{A n B). Moreover, letting 
/i denote the Lebesgue measure on , we define the symmetric difference distance, also called the symmetric 
distance for simplicity, c?0 : C x C — > K>o by 

de{A,B) = ^i{AQB), 
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that is, the symmetric distance between two sets is the measure of their symmetric difference. Given these 
notions, it is useful to identify sets that differ by a set of measure zero: we write A ^ B whenever fi{AQB) = 0. 
Clearly, ^ is an equivalence relationship on C and, accordingly, we let C* ~ C/ ^ denote the quotient set of 
dosed subsets of Q. Now, for any two elements A* and B* of C*, we define dQ{A*,B*) = dQ{A, B) where A 
and B arc any representatives of A* and B* , respectively. With this notion of dQ on C* xC*, it is easy to verify 
that (C*,(ie) is a metric space. However, to the best of our knowledge no compactness result is available for 
(C*,rfe). 

Next, we introduce a particular family of subsets of C whose structure is sufficiently rich for our algorithm. 
For ^' G N, let C(^) C C denote the set of i-convex and closed subsets of Q, that is, the set of subsets of Q equal 
to the union oi l convex and closed subsets of Q. Formally, we set 

C(^) = {w C Q I u = U^^]^ Si where Si, . . . , Si C Q are convex and closed}. 

Note that we do not require the sets Si, . . . , Se to he distinct so that C(fc) C C(|), for any k < i. In what follows 
we study the quotient set of i-convex and closed subsets C^*^j = Cj^)/ ~. The next result is the main result of 
this section. 

Theorem 5.1 (Compactness of C^*^^) The pair (C*^^,c?e) is a metric space and, with the topology induced by 
dQ , the set C*^-^ is compact. 

Proof. It is easy to verify that is a metric on C* . Instead, proving the compactness of Cj*^j requires some 
attention. We aim to show that any sequence in C*^^ has a subsequence that converges to a point in C^^^. This 
fact's proof is articulated in several steps and relies upon several known results: 

(i) the space C, endowed with the Hausdorff distance dn : C x C ^ R>Oj is [301 Theorem 0.8] a compact 
metric space; 

(ii) if a sequence of closed convex subsets of Q converges in the Hausdorff sense to a set K, then [531 Propo- 
sition 1.6.8] K is closed and convex; and 

(iii) for any two convex subsets A, _B of Q C M'', it is known [TSl Eq. (1)] that 

where D = max{diam(A), diam(i?)} and where is the volume of the unit sphere in W^. 

Now let {u*(n)}„gz>(, be any sequence in C^*^^. For all n G Z>o, pick any representative of the equivalence class 
v*{n), denote it by v{n) and consider the sequence {w(n)}„gz>Q. Because {w(n)}„gz>o is a sequence in C(^) and 
because C(£) C C, it follows from fact (i) above that {w(n)}„gz^p contains a subsequence that converges in the 
Hausdorff sense to an element of C. In other words, there exist {'y('T-fc)}fcez>o ^ {''^('^)}nez>o and w G C such 
that 

lim v{nk) = V, 

k — >oc 

where '= denotes convergence in the Hausdorff sense. We claim that v G C(f) so that the set C*^^ is compact. To 
show this claim, we plan to exhibit a collection of convex and closed subsets of Q, say {^i, . . . , Si\, such that 

We begin to prove this claim as follows. By definition of C(f), for all k G Z>o there exists a collection 
{5i(nfc), . . . , 5'f (nfc)}, of convex and closed subsets of Q whose union equals v{n}^). Now, we consider the 
sequence {Si{nk)}k^z>o- Again, since (C, d^j) is a compact metric space we have that there exists a subsequence 

above ensure that Si is a 

convex closed subset of Q. Consider now the sequence {>S'2(fifci)}fciez>o. By reasoning as previously wc know 
that there exists a subsequence {"■fc2}fc2GZ>o ^ {"'fei}fciGZ>o such that lim^^^oo S2{nk2) — ' S2 where 5*2 is some 
convex closed set of Q. Moreover, it is clear that also limfc2_oo 5*1 (nfc^) '= 6*1. By iterating this procedure we 
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{"-fci}fciGZ>o ^ {nk}kei.>o such that linik^^oo Siiuk^) = S'l for some ^i. Fact 



conclude that there exist a sequence {ns}sez>o and a collection of convex and closed sets {^i, . . . , Si} such that 

linis^oo Si{Tis) = Si for aU i G {1, . . . ,£}. 
Next, we aim to show that 

e e e 

lim I J S,{n,) = 11 lim S,{n,) = I J 5,. (15) 

Z=l 2=1 i—1 

For simplicity, let us denote lims_too Ui=i *S'^(f^s) by S'oo- For p ^ Q, note 

p e 5oo dist ^p, ^hm^ [J 5*^(71,5 )^ =0 <^==> ^liiii^ dist ^p, [J 5, (713)^ =0, 

where, for a given closed set X, dist(p, X) denotes the Euclidean distance between p and X, namely, minajgjf \\p— 
x\\. Using the fact that, for given closed sets X and Y, dist(p,XUF) = min {dist(p, X), dist(p, F)}, we can 
write 



hm dist ^p, IJ 5,(71^)^ = 



lim min {dist(p, Si{ns)), . . . , dist(p, Si{ns))} = 

s — >oo 

3j e {!,...,£} s.t. lim dist(p, 5,(ns)) =0 

3je {!,...,£} s.t. peSj 

e 



The above chain of implications proves pS]) . Now observe that, from the uniqueness of the limit it follows that 

limfc^oo w('T-fc) = lims^oo Ui=i "^il^-s) ~ Uj=i "^i ~ ^- Since Si is closed and convex for all i G {!,...,£}, it 
follows that V G C(^) and, in turn, that C(^) endowed with the Hausdorff metric is a metric compact space. 

To establish the statement of the Theorem it finally remains to prove that either limfc__»oo '^e(^("-fc), v) = 
or lims^oo dQ{v{ns),v) = 0. To this end, observe that, given Xi, X2,Yi,Y2 C Q, the following inclusion holds 
{Xi U X2) e (Yi U Y2) C {Xi e Yi) U {X2 e Y2). Hence, we compute 

V (n,) ev=([j S,in,) \ Qv=([j S^ns)] Q (\J s}j C (J e 5.) , 

\i=l / \i=l / \i=l / 1=1 

which implies 

de (w(ns) ,t;) = de I U Si{ns),v\ <^dQ {S,{ns),Si) 

\i=l ) i=l 

J-^ 2Kd /max{diam(S'^(7is)) ,diam(S'i)}\ , 
^ oi/rf _ 1 ^ ^« (^^("«)' 



- 2i/rf - 1 \^ 



where the second and the third inequalities follow, respectively, from fact (iii) above and from the upper bounds 
max{diam(5i(ns)) ,diam (^i)} < diam((5), for all i G {1,...,^}. Since liuis-too dn {Si{ns), Si) ~ for all 
i G { 1 ,...,£} , we conclude that 

lim dQ{v{ns),v) = 0. 

s — *oo 

Now, let V* denote the equivalence class for which w is a representative. Since the metric rfe is insensitive to sets 
of measure zero, it follows that Imis^oo dQ(v* {ns),v*) ~ and, in turn, that {?;*(n)}„gz>o has a subsequence 
convergent to point of C^iy This concludes the proof that C*^^ is a compact space. I 
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The metric cIq naturally extends to a metric over the space (C*)^, and hence over (C?^n)^, by defining 



N 



dQ{u,v) ^^dQ{ui,Vi), (16) 
1=1 

for any u = and v = ('fi)i^i hi (C*)^. The compactness of (Cf^J^ is then a simple consequence of 



Theorem 15.1 



Corollary 5.2 (Compactness of (C^*^))^) T/ie pair ((C*^))^, cJq) is a metric space and, with the topology 



induced by dQ, (C,*^-,)^ is compact. 



6 On the continuity of some geometric maps and the resulting con- 
vergence proofs 

In this section we prove the main convergence theorems for our gossip coverage algorithms. First, however, we 
need to establish the continuity properties of certain geometric maps and of the proposed algorithms T and . 
Before proceeding, we discuss two significant modeling aspects. First, recall that SectionOintroduces the spaces 
C, C(£), C^, C^-^ and the corresponding quotient sets C*, C*^y (C*)^, (C*^^^. We can do the same with the 
partition space Vn ■ Indeed from Defiiiitionl2.]tiii) we have that each component tj^ of u G Vn can be mapped by 



the canonical projection into an equivalence class v* in C* \ {0}; in turn, any v S Vat can be naturally associated 
to the iV-collection of equivalence classes v* = (w*)fli . Accordingly, we denote the space of equivalence classes 
of N -partitions by Vj^. Recall that all these quotient spaces are metric spaces with respect to the symmetric 
distance dQ. 

Second, recall the following maps: the centroid map Cd : {A £ C \ n{A) > 0} ^ Q defined in equation ([5]), 
the 1-center function Tii : Q x C ^ K>o defined in equation the multicenter function Ticentroid '■ Vat —^ IR>o 
defined in equation p^ . the gossip coverage map Tij : Vn — > Vat with i ^ j defined in equation ([TT|) . and, 
for (5 > and i ^ j, the modified gossip coverage map Tf^ : Vn — > Vn defined in Section [3^ We claim that 
all these maps are insensitive to sets of measure zero. To substantiate this claim, observe that the integrals 
of a bounded measurable function over a set A and over a set B are equal if dQ{A,B) = 0. This observation 
allows us to redefine the centroid map, the 1-center function and the multicenter function as Cd : C* \ {0} — > Q, 
Til : Q X C* ^ IR>o and Hccntroid : — > M>o, respectively. Regarding the modified gossip coverage map Tfj, 
we reason as follows. For v* G V^, let v G Vn and v' e V^v be two representatives of v* and let v and v' denote, 
respectively, the images of v and v' under the map Tfj, that is, v = Tfj{v) and v' = T-j{v'). Since the centroid 
map and the definitions of the points pi and pj in equation (|14p are insensitive to sets of measure zero, it follows 
that c?e(w, u') = 0; in other words v and v' belong to the same equivalence class, say v* G V^. From these facts 
we can redefine the modified gossip coverage map as : — > V^. An analogous argument applies to the 
map T. This concludes the justification of our claim. Finally, note that the Voronoi map V : \ Sn — »■ Vn 
defined in equation ^ can be composed with the standard quotient projection map and therefore denoted by 
V : Q^\Sn ^ Vj^. In summary, the centroid map, the 1-center function, the multicenter function, the modified 
gossip coverage map, and the Voronoi map are indeed insensitive to sets of measure zero. 



6.1 Continuity of various geometric maps 

We start by recalling that the compact connected set Q is equipped with a bounded measurable positive function 
(j> : Q ^ K>o- We define the diameter of Q and the infinity norm of 4> by diam((5) = max{||a; ~ y\\ \ x,y & Q} 
and ||c/'||oo = ^SiXx^Q 4>{x), respectively. The following lemma states some important properties of the 1-ccnter 
cost function. 

Lemma 6.1 (Continuity properties of the 1-center function) Given a compact convex set Q C M'*, let 
(j) : Q M>o be bounded and measurable and let f : ]]i>o — > IR.>o he locally Lipschitz, increasing, and convex. 
Define the function Hi : Q x C* ^ ]R>o as in equation ([?]). The following statements hold: 

(i) the function p ^ 7ii(p; A) is strictly convex in p, for any A £ C* \ {%}, 
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(ii) the function p t-^ 'Hi{p; A) is globally Lipschitz in p, for any A G C* , and 
(Hi) the function A 'Hi{p; A) is globally Lipschitz in A with respect to eIq, for any p £ Q. 
We now state the main result of this subsection. 

Theorem 6.2 (Continuity properties of centroid, multicenter and Voronoi maps) Given a compact con- 
vex set Q C , let (f> : Q ]R>o be bounded and measurable and let f : M>o — > K>o be locally Lipschitz, 
increasing, and convex. With respect to the topology induced by d^, , the following maps are continuous: 

(i) the centroid map Cd : C* \ {0} —t Q, 

(ii) the multicenter function Ticcntroid '■ '^n ^ ^>o, 

(Hi) the Voronoi map V : \ Sn Vj^, 

(iv) for all i,j £ {1, . . . , N}, i ^ j , the gossip coverage map Tij : {v G V'^ \ Cd{vi) ^ Cd{vj)} — > V^, and 



(v) for all S > 0, i,j£ {1, . . . , N}, i ^ j, the modified gossip coverage map Tfj : 



'AT- 



The continuity properties (ii) and (iv) (respectively, (v) ) are exactly what is needed to apply the Krasovskii- 
LaSalle invariancc principles stated in Section 3] to the gossip coverage algorithm (respectively, to the modified 



gossip coverage algorithm). The continuity properties (i) and (iii) are intermediate results of independent 
interest. 

Proof. [Proof of Lemma [^T] Let Lf he the Lipschitz constant of / : [0, diam(Q)] — > R>o. We check the claims 
in order, beginning with statement (i) For A G (0, 1), A G C* \ {0}, and pi,p2 S Q, we compute 



HiiXpi + {1 - X)P2; A) 



fi\\Xp, + {l-X)p2-qmq)dq 



< / f{X\\pi^q\\ + il-X)\\p2-qmq)dq 
JA 

< 



(17) 
(18) 



(a/(||pi - gil) + (1 - A)/(||p2 - q\\))<f>iq)dq 
= A7^i(pi;v4) + (l-A)Hi(?52;A), 

where inequality (jl7[) follows from the triangle inequality and from / being increasing, and inequality (|18p 
follows from the convexity of /. This inequality proves convexity. Moreover, since the first inequality is strict 
outside the line passing through pi and p2 and since A has non-empty interior, the function is in fact strictly 
convex. Note that statement (i) implies that p i— > T-Ci{p;A) is locally Lipschitz, using [3H Theorem 10.4]. The 



stronger statement (ii) can be derived as follows. For pi,p2 G Q, we compute 



\nApi;A)-Hi{p2;A)\ 





1 

JA 


- q\\)(f){q)dq - 


[ f{\\P2~qmq)dq 
'A 




[ IfiWpi 

JA 


-'7ll)-/(lb2- 


-qm{q)dq 






f l/(lbi 

A 


-'zll)-/(lb2- 


-qm{q)dq 






f %lbi 

A 


-P2\\<P{q)dq < 


LfUWoc^KMPl ~ P2 
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This implies the Lipschitz condition in statement (ii) Statement (iii) can be proved as follows. Let A, A' be 
two elements of C* , note A {A\ A')[j{An A') and compute 



\niip;A)-niip;A')\ = 



A\A' 



.f{\\p-qmq)dq 



A'\A 



f{\\p-qmq)dq 



< 



A\A' 



.f{\\p-qmq)dq 



f{\\p-qmq)dq 



f{\\p- q\\)<P{q)dq - 

A\A' J A'\A 



A'\A 

f{\\p^q\\)mdq 



f{\\p-q\\)mdq 



AQA' 



< max{/(|[p -q\\)\p,qe AqA'} ||0||oo 6 A') 

< /(diam(g)) II0II {A, A'), 

where last inequality follows from / being increasing. The bound implies the Lipschitz condition. 



Before proving Theorem 16. 21 we need the following lemma about perturbations of convex optimization prob- 
lems. 

Lemma 6.3 Given a compact convex set Q C M.'^ and a metric space {X,d), let H : Q x X ^ M. have the 
properties that 

(i) the map x > H{q,x) is globally Lipschitz for all q £ Q, and 
(ii) the map q ^ H{q,x) is continuous and strictly convex. 
Then the map q* : X Q, defined by q*{x) ~ argmiUggg _ff (q, a;) , is continuous. 

Proof. Let be the Lipschitz constant of x ^ H{q,x). Thanks to the Lipschitz condition of the function 
X y-^ Hiq, x), for all x, y £ X, the point of minimum q*{x) takes value in S = {q G Q \ H{q, y) < H{q*{y),y) + 
2L//(i(?/, x)}. Since 5* is a sub-level set of the strictly convex function q i— > H{q,y), and since the diameter of 
a sub- level set depends continuously on the level, the distance — can be made arbitrary small by 

reducing d{y,x). This implies the claimed continuity. I 

We are now ready to prove the main result. Proof. [Proof of Theorem 16. 2j We prove the theorem claims in 

the order in which they are presented. Claim (i) follows combining Lemma 16.31 and Lemma 16.11 Since the 
multicenter function Ticcntroid is a sum of suitable 1-center functions "Hi, the claim (ii) is also immediate. 



Regarding claim (iii) we discuss in detail the two dimensional case. We first prove the continuity when 
N = 2: let pi and p2 be points in Q. Let 21 = \\pi — P2\\- Since pi ^ p2, I > 0. Up to isometrics, we can 
assume that, in the Euclidean plane {x,y), pi = {—1,0) and p2 = {1,0)- Let di and d2 be the distances from 
the origin of points pi and p2, respectively. It is clear that the two Voronoi regions of pi and p2 are separated 
by the locus of points {x € Q \ \\x — pi\\ = \\x — P2\\}, that is the vertical axis. Now, we assume that the 
positions of pi and p2 are perturbed by a quantity less than or equal to e, with < e < L By effect of the 
perturbation, the axis separating the two Voronoi regions is perturbed, but it is contained in the locus of points 
^12 (e) = {x Q\ ||x — Pill — \\x — P2II ^ 2e}. By definition, this is the set comprised between the two branches 
of the hyperbola whose equation is — 73377 ~ 1. By elementary geometric considerations, the area of this 
region can be upper bounded by 



M(yi2(e)) < 2e2 diam(g) + 4diam(Q)^ 



e/l 



< 4 diam(g)(l 



1 - ^ 



diam((5) 
I 



This bound implies the continuity. The case in which N > 2 follows because, moving all points by at most e, 



the change in all the regions is upper bounded by IJj^ 



<i.j<N -'y 



Yij{e), which vanishes as e — » 0+. 
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The last remaining step is to prove claims (iv) and (v) We focus on claim (v) and show claim (iv) as a 
byproduct. Let v G Vjv and let v = Tfj{v). According to the definitions in Subsection 13.41 v is characterized by 

the sets Ri and Rj . Recall that these two sets depend on the sets Ri Rj , on the scalar /3ij {v) and on the points 
Pi and pj. One can see that Pij{v) is a continuous function of its arguments Vi and Vj. Hence, it suffices to 
show that also Ri, Rj and pi, pj depend continuously on Vi and Vj. To do this, introduce v' £ Vn and compute 
v' — Tfj{v'). Assume Gd{vi) ^ Cd(uj) and Gd{v[) ^ Cd(Uj). Analogously to how we defined Ri and Rj, we now 
define the regions R'^=v[r\ i7bisoctor(Cd(?;^), Cd(w^)) and R'^ = v'^ n i/biscctoi(Cd(u-), Cd(u^)). We aim to upper 
bound the composite distance dQ{Ri,R[) + dQ{Rj, R'^). Observe that this composite distance depends on the 
difference of the two argument regions and V j , v'- both directly and indirectly via the induced difference 



between the centroids. Recalling the proof of claim (iii) let e be an upper bound on the displacement between 
the two centroids. Then the region Y{e) = {x e Q | \\\x — Gd{vi)\\ ~ |jx — Cd(wj)|| < 2e} needs to be included 
in the upper bound on the composite distance. Combining these considerations we obtain 

dQ{R„R!,) + dQ{R„R!j) < {de{v„v',) + de{vj,v'^)) 

+ M(F(max{|l CdK) - Cd{v[)\\, \\ Cd{v,) - Cd(«;)||})). (19) 

Clearly, if deiv„v^) and de{vj,v'^) -> 0, then || Cd(i;,) - Gd{v^)\\ ^ and || Cd(v,) - Cd(w,')|| ^ and, in 
turn, also dQ{Ri,R[) + dQ{Rj, R'j) — ^ 0. Hence, we can argue that the sets Ri and Rj depend continuously on 



the regions Vi and Vj . This is enough to prove statement (iv) provided the two regions Vi and Vj have distinct 
centroids. Moreover a direct consequence of this fact is that also the points pi and pj depend continuously on 
Vi and Vj. Finally, observe that in the limit case Gd{vi) = Cd(z;j) the continuity of Ri and Rj is captured by 
the fact that (iij[v) is a continuous function of Gd{vi),Cd{vj) and that (3ij{v) = if Gd{vi) = Gd{vj). The 
continuity of Ri, Rj, I3ij{v), pi and Pj imply also the continuity of Ri and Rj and, in turn, of T^-. I 



6.2 Convergence proofs 



In view of the identification between A^-partitions and their equivalence classes introduced at the beginning of 
this section, we are now ready to complete the proof of the convergence results presented in Section [3.31 

We start by clarifying the precise meaning of convergence in Theorems 13. 71 and 13.81 Specifically, we say that 
a sequence of partitions {w(0}tez>o C V^v converges to a set of partitions X C Vn if the symmetric distance 
from {v{t)}te'L>o *° ^ converges to zero, that is, 

lim inf{de(t;(0, x) | x G X} = 0. 

t — >oo 

Proof. [Proof of Theorem 13. 8j We prove the deterministic statement (i). We start by observing that, through 
the canonical projection, the evolution {i'(i)}t6Z>o C Vat of T* : Vjv Vn can be mapped into the evolution 

{^^*(i)}tGZ>o C of T'' : ^ V]^. We aim to" apply Theorem to the dynamical system 



' : V 


N ~^ 




(i) 


(ii) 




(iii) 



and its evolution {w*(t)}tez 
and (iv) of Theorem 14.31 are satisfied 

Since {v(t)}ti£z,>o non-vanishing and finitely convex by assumption, it follows that there exists £ G N such 
that the fi-limit set of {'i'(i)}t6Z>o is contained in C(£) fl Vn, that is, Q{v{t)) C C(f) n Vat. This imphes also 
that i}{v*{t)) C C(*^^ n Vn- As stated in Theorem 15.11 C*^^ is compact in the topology induced by the metric 
Hence, even though C*^-, is not strongly positive invariant for , the weaker version of Assumption (i) of 



Now, as one can deduce 



"e- even tllUUgll l^^^j 

Theorem I4.3[ as given in Remark 14.61 holds true for the sequence {v*(t)}tgz> 

from Theorem I6.^(ii)| and Lemma 13.41 the function Ticcntroid satisfies the Assumption (ii) of Theorem | 
thus playing the role of a Lyapunov function for the dynamical system . Moreover, from Theorem I6.^(v)[ 
note that the system evolves through maps that are continuous in with respect to the metric de- thus 

is satisfied. 



the Assumption (iii) of Theorem 14.31 is satisfied. Finally observe that the Assumption (iv) of Theorem 14.31 
corresponds to the assumption of uniform persistency in Theorem l3.81 Therefore, we conclude that the evolution 
{v*(i)}*6Z>o converges to the intersection of the fixed points of the maps Tfj, for all i,j G {1, . . . ,N}, j ^ i. 
According to Lemma 13.31 this intersection coincides with the set of mixed centroidal Voronoi partitions up to 
sets of measure zero. 
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The proof of the stochastic statement (ii) foUows the same hnes, applying Theorem 14.51 instead of Theo- 
rem US] ■ 



Proof. [Proof of Theorem 13. 7] The proof of this rcsuh follows the lines of the proof above with two distinctions. 
The distinctions come from the assumption that the evolution v, in addition to being non- vanishing and finitely 
convex, is also distinct centroidal. The first distinction is as follows. In order to apply the Krasovskii-LaSalle 
invariance principle we require the continuity property stated in Theorem I6.1^(iv) Additionally, we note that 
the space of finitely convex and distinct centroidal partitions 

{v e Vn nCjJ) I II Cd{v,) - Cd{vj)\\ > e for all i j} 

is a closed and hence compact subspacc of Vn- The second distinction is as follows. Since we rule out the case of 
coincident centroids, we can infer convergence to centroidal Voronoi partitions instead of convergence to mixed 
centroidal Voronoi partitions; see Lemma 13. t^(ii)[ I 



7 Conclusions 

In summary, we have introduced novel coverage, deployment and partitioning algorithms for robotic networks 
with minimal communication requirements. To analyze our proposed algorithms, we have developed and char- 
acterized (1) intuitive versions of the Krasovskii-LaSalle Invariance Principle for deterministic and stochastic 
switching systems, (2) relevant topological properties of the space of partitions, and (3) useful continuity prop- 
erties of a number of geometric and multicenter functions. 

We believe there remain interesting open issues in the study of gossiping robots and of dynamical systems on 
the space of partitions. We are keen on extending these ideas to non-convex complex environments and discrete 
environments such as graphs. Following Remark 13.21 we plan to study gossip coverage algorithms for more 
general multicenter functions, including nonsmooth, anisotropic and inhomogeneous functions. Additionally, we 
plan to investigate gossip coverage algorithms capable of adapting to time- varying scenarios such as problems in 
which robotic agents arrive to and depart from the network. Finally, inspired by stigmergy in territorial animals, 
we plan to design communication protocols for multiagent systems based on the ability of leaving messages in 
the environment. 



A A robotic network implementation of gossip coverage algorithms 

In Subsection 12.21 we discussed coverage control algorithm for groups of robots with synchronized and reliable 
communication along all edges of a Delaunay graph; then in Subsection 13.21 we introduced our gossip coverage 
algorithms. Here we discuss in full detail one possible way of implementing the partition-based gossip coverage 
algorithm in a robotic network with weak communication requirements. 

We consider a group of agents all having the following capabilities: (CI) each agent i G {1, . . . , N} knows its 
position and moves at positive speed Ui to any position in the compact convex environment Q C M^; (C2) each 
agent may store an arbitrary number of locations in Q and has a clock that is not necessarily synchronized with 
other agents' clocks (specifically, we assume same clock skew, but different clock offsets among the agents); and 
(C3) if any two agents are within distance rcomm of each other for some positive duration of time, then they 
exchange information at the sample times of a Poisson process with intensity Acomm- 

The random destination & wait algorithm is described as follows. Given a parameter e < rcomm/4, each agent 
z S {1, . . . , A^} maintains in memory a dominance region Vi and determines its motion by repeatedly performing 
the following three actions over periods of time that we label epochs: 
1: instantaneously agent i selects uniformly randomly a destination point qi in the set {(/ £ | dist {q,dvi \ 

dQ) < e}; 

2; agent i moves in such a way as to reach point qi in time precisely equal to d = diam((5)/ minjwi, . . . , u^}; 
and 

3: agent i waits at point qi for a duration that equals either d with probability 1/2 or 2d with probability 1/2. 
We have assumed that agents may reach locations at a distance up to e away from Q; this assumption can be 
removed at the cost of additional notation. 
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The random destination & wait algorithm is meant to be implemented concurrently with the modified 
gossip coverage algorithm - where the parameter 6 is selected to satisfy 5 < rcomm/4. The two algorithms 
jointly determine the evolution of the agents positions and dominance regions as follows. If any two agents 
i and j are within communication range rcomm at any instant of time and for some positive duration, then, 
at each sample time of the corresponding communication Poisson process, the two agents exchange sufficient 
information to update their respective regions Vi and Vj via the modified gossip coverage map Tf^ . 

Proposition A.l (Random destination & wait ensures persistent random gossip) Consider a group 
of agents with capacities (CI), (C2) and (C3) and parameters Ui, rcomm o-nd Acomm- Assume the agents im- 
plement the random destination & wait algorithm and the modified gossip coverage algorithm with parameter 
e < ''comm/4 and 5 < rcomm/4. The following statements hold: 

(i) the sequence of applications of the modified gossip coverage map is a randomly persistent stochastic process; 
and 

(it) the resulting evolution v : R>o — *■ Vat, conditioned upon being non-vanishing and finitely convex, converges 
almost surely to the set of mixed centroidal Voronoi partitions. 



Proof. Assume that at time s £ K>o, the regions Vi and Vj are at distance S or less. Pick any two points 
Pi € dvi \ dQ and pj S dvj \ dQ such that \\pi — Pj\\ < 2(5 and define two balls of radius e centered at pi and pj 
by Bi = {q & R"^ \ \\q — Pi\\ < e} and Bj = {g g | \\q — pj\\ < e}. By construction, each point in Bi is at most 
at distance 2e + 2S < rcomm from each point in Bj. We claim that, independently of the agent positions and 
the environment partitions at time s and of their evolution throughout the past time interval [0, s], there exists 
a duration A £ K>o such that the meeting event "agent i is in Bi and agent j is in Bj during the time interval 
[s + A, s + A + c?]" has probability that is lower bounded by a positive constant. If these meeting events have 
positive probability, then communication events have positive probability and therefore the proposition follows. 

We prove the claim as follows. At time s, we say that an agent is in OFF state if it is at the initial time of an 
epoch or it is moving to its next destination point. We say it is in ONi state if it has just reached a destination 
point or it has been waiting at the destination point for a duration of time less than d. We say it is in ON2 
state if it has been waiting at the destination point for a duration of time greater than or equal to d and less 
than 2d. Note that each agent remains in each one of the three states {OFF, ONi, ON2} for a duration of time 
d. The transitions between states are regulated by a Markov chain with three states {OFF, ONi, ON2} and with 
(column-stochastic) transition matrix 

"0 1/2 1' 
A= 1 . 
_0 1/2 0_ 

For now, suppose that at time s, agent i is at the initial time of an interval of duration d in either state OFF, or 
state ONi or state ON2. From the matrices 



A^ 



1/2 
1/2 




1/4 
1/2 
1/4 



1/2" 

1/2 



and A^ 



1/2 3/8 1/4" 
1/4 1/2 1/2 
1/4 1/8 1/4 



we establish the following three facts: 

(i) if agent i is at state OFF at time s, then at time s + 3d it is at state OFF with probability 1/2 and at time 
s + 5(i it is at state ON2 with probability 1/4; 



(ii) if agent i is at state ONi at time s, then at time s + 3d it is at state OFF with probability 1/4 and at time 
s + 5d it is at state ON2 with probability 1/8; and 

(iii) if agent i is at state DN2 at time s, then at time s + 3d it is at state OFF with probability 1/2 and at time 
s + 5d it is at state ON2 with probability 1/4. 

We now return to the two agents i and j. According to the robot capability (C2), at time s the clocks of 
agent i and agent j, denoted by ti{s) and tj{s) respectively, are distinct and can be written as an integer 
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number of durations d plus a quotient; specifically ti(s) = kid + rji, and tj{s) = kjd + rjj, with ki,kj G Z>o, 
771,772 G [0,(i[, and possibly /c,; 7^ kj and 77^ 7^ 77^. Without loss of generality, assume 77,; < 77^. Facts (i), (ii), 
and (iii) and the assumption that clock skew is equal among agents, jointly imply that, with probability greater 
or equal to (1/8)^ = 1/64, agents i and j during the interval [s + 4(i — 77^, s + 5(i — 77^] are waiting at some 
destination point chosen, respectively, at time s + 3d ~ rji and s + 3d — r]j. If agent i is in Bi and agent j is 
inside Bj during the interval [s + 4d — 77^, s + 5c? — 7;^], then they communicate at least once with probability 
Pcomm := 1 — e~'^'^'=°'"" > 0. Let us now find a lower bound on the probability of agents i being in Bi at time 
s + Ad — rji] clearly this happens if the following two facts occur: 

• during the interval [s, s + 3d — rji] agent i does not communicate with any other agent and hence its 
dominance region Vi remains unchanged; and 

• at time instant s + 3d — rji agent i is at the beginning of a new epoch (state OFF) and selects a destination 
point in Bi. 

Since during the interval [s ^ rji, s + 3d — rji] agent i spends a time duration of at most 2d at distance less 
than 6 from any other agents, then the probability that agent i does not communicate in that interval of 
time is lower bounded by (1 — Pcomm)^ = e"^'^'''"'""". Moreover, provided that at time instant s + 3d — i^i 
agent i is at the beginning of a new epoch (state OFF) and that the Vi is unchanged until time instant s, 
we have that the destination point for agent i is selected uniformly randomly so that it belongs to Bi with 
probability area(i3i)/ area ({g £ | dist [q,dvi \ dQ) < e}) > Parca 7re^/( area((5) + e pcrimctcr(Q)) . 
Similar considerations hold also for agent j. 

In summary, if agents i and j are neighbors at time s, then they communicate at least once during the interval 
[s + 4d — rji, s -\-bd — rji] with probability lower bounded by (l/64)p^^(,jj(l — Pcomm)^Pcomm > 0. If agents i and 
j are at distance greater than 5 at time s (and therefore they are not neighbors), then the modified gossip 
coverage map is the identity map so that their communication is immaterial. I 



B A counterexample showing the necessity of uniformly persistent 
switches 
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Theorem 14.31 contains a persistent switching conditions, that is, it requires the existence of Z) G N such that 
every map T^, i G {1, . . . , 771}, is applied at least once within every interval [n, n + Z?[, n G Z>o. This appendix 
contains an example proving the necessity of this condition. 

Consider the plane in polar coordinates X = IR>o x [0, 27r[U{0, 0}. Define the standard metric d : X x X 
M>o as follows: let (pi,0i), (^2,^*2) be any pair of elements of X and 

d((pi,6'i), {p2,d2)) = \/(/9i cos 6*1 - p2 cos 6*2 ) ^ + (pi sin 6*1 - p2 sin 6*2 ) 
Consider now the continuous maps Ti : X ^ X , i £ {1, 2}, defined by respectively 

Ti{p,0) = 




p — 1) mod 27rj , if p > 1, 

((1 - sin 6i)p, 61), ifO<6l<7r, 
(p,6l) if7r<6i<27r. 



Define T : X ^ X hy T{p, 9) = {Ti (p, 6) , T2 (p, 9)} and the function U : X ^ R>o hy U (p, 6) = p. Observe 
that U is continuous and non-increasing along T. Assume now that there exists Z? G N such that, for any 
n G Z>o, there exist ni and 772 within the interval ]n,n + D] such that Xm+i = Ti{xni) and Xn2+i ~ T2{xn2)- 
Then, by Theorem 14.31 the w-limit set of each evolution of T is a subset of 

{(p,6l) G X| p = l,7r < 61 < 27r}U {0,0}. (20) 

Next, we relax the condition that the map T2 is applied at least once inside each interval of arbitrary ampli- 
tude D and we show that there exists one sequence that does not converge to the w-limit set in equation (|20p . 
To this aim, assume the sequence {{p{n),6(n))}n^z>a satisfies 
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(i) p{0) > 1; 



(ii) (p(l),0(l)) = Ti(p(O),0(O)) and 

(iii) {p{n+l),e{n + l)) = Ta {p{n),e{n)) if and only if tt < 9{n) < 27r and {p{n),e{n)) = Ti - 1), 6'(n - 1)). 

Note that if tt < 9{n) < 2tt, then T2 {p{n),9{n)) = {p{n),6{n)). Therefore, the evolution {{p{n),9{n))} equals 
|(p(n),^(ri))| where (p(0),?(0)) = {p{0),e{0)) and T^" (p(0), ^(0)) . Regarding this new se- 

quence, observe that 

1 < p{i) < 2 and p(i + 1) < p{i), for all i > 1, (21) 
0<9{i + l)- 9{i) < TT, for aU i > 1, and lim (9{i + 1) - 9{i)) = 0, (22) 

r r r y \ ~^ 

}}^ E (^(* + 1) - ^(*)) = E - 1) = }}^ E ^YTY^ + ' " = ^' ^^^^ 

where the equality p{i) — 1 = (^^^j^ + i — 1)^^ can be proved by induction over i. Properties ([2T|) . ([22|) . 
and ([23]) ensure that there exists a sequence {n^ | /i G Z>o} such that (p{nh), 9{nh)) = T2 {p{nh — 1), 9{nh — 1))) 
for aU h € Z>o, and {p{t),9{t)) = Ti {p{t - l),0(t - 1))) ii t ^ {nh \ h G Z>o}. Moreover, we have that 
lim;i_oo {nh+i—nii) = 00. In other words, both the maps Ti and T2 are applied infinitely often along the evolution 
described by {{p{n),9{n))}, but there does not exists D €N such that T2 is applied at least once within each 
interval [n, n+D], n G Z>o. Observe that, in this case, {p{n), 9{n)) converges to the set {(1, 9)\9 G [0, 2it]} C X. 
This set is different from the cj-limit set in equation ([20|) . 



C Discontinuity of the multicenter function in the Hausdorff metric 

To see that Ticontroid '■ V]^ M>o is not Hausdorff-continuous, consider the sequence of 2-partitions {f (t)}tGZ>o 
of the interval [—1,1] CM defined by 



-1,-1 



2*+i 



2'-^-l 

u u 

ft = -(2*-i-l) 



1 h 



2t~i 2*+!' 2*-i 2*+i 



U 



1 - 



9t+l ■ 



and by V2{t) = [—1, 1] \ vi{t). Note that both sequences {vi{t)}t^z>a and {v2{t)}tez>o converge to [—1, 1], and 
that Cd{vi{t)) = Cd(w2(0) = Cd([-1, 1]) = 0, for all t > 0. Hence,~for (j){s) = 1 and /(s) = s, we compute 



niiO,viit)) = [ \x\dx^2 j I 
Jvi(t) Ji,i(t) n[o,i] 



xldx = 2- = 1, 
' 2 



and consequently Hcontroid('^(0) — 2, while 7icontroid(liint^oo i^('^)) = 27-^i(0, [— 1, 1]) = 4. This shows that 

limt^oo ■Wccntroid(w(t)) 7^ 'Hccntroid (hmj^oo w(i ) ) • 
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